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Statistical mechanics of spin models with disorder is an important field of condensed 
matter theory onto which much attention is focused at present time. All experi- 
mentally accessible magnetic systems with localized magnetic moments are always 
affected by randomness and the understanding of disorder effects within the frames 
of simple models can be of much use in analysis of experimental data. Usually two 
types of disorder are considered, namely, the diagonal and off-diagonal disorder (see 
for example [1]). However, it is evident that for alloys as well as for topologically dis- 
ordered systems the random diagonal and off-diagonal matrix elements involved into 
the Hamiltonian of the system are connected and one faces the correlated diagonal 
and off-diagonal disorder. Starting from 1974 there exists a model of tight-binding 
electrons with correlated Lorentzian disorder the density of states for which can be 
found exactly [2] . Latter on the method elaborated in Ref . [2] was used for obtaining 
the density of magnon states for the spin-i XX chain with correlated Lorentzian 
disorder [3]. However, the obtained analytic results pertain only to thermodynamics 
and an interesting question about spin correlations and their dynamics still remains 
open. The aim of the present paper is to discuss this problem using a numerical 
approach developed recently [4,5]. 

We consider a chain of spins | governed by the Hamiltonian 
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The exchange couplings are assumed to be random independent variables each with 
the Lorentzian probability distribution 
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where Jq is the mean value of exchange coupling and F controls the strength of 
disorder. The transverse fields are assumed to be related to the exchange couplings 
in the following way 
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where Qq is the mean value of transverse field at site. We are interested in the zz 
time-dependent correlation functions (s|(t)s|_,_^) and the random-averaged trans- 
verse dynamic structure factor and susceptibility defined by the formulae 



n=l 



N 



" ■ fOO ■ ■ — p = — 

Jo I U J i 

respectively. 

An idea of numerical computations of these quantities is explained in details 

in Rcfs. [4,5]. We considered chains of = 300 spins with Jq = —1, Qq = 0.5 
and r = 0.1 at the low temperature j3 = 1000. We computed correlation functions 
{sl5o{t)450+n) - (sf5o)(si50+n) with n = 0, ±1, . . . , ±100 for the times up to ^ = 200, 
put e = 0.01 and averaged the zz dynamic structure factor and susceptibility over 
3000 random realizations. To understand the accuracy of the obtained results we 
performed additional calculations revealing the effects of finite N and n for the 
considered times t at the taken value of temperature p. We also compared the 
numerical results with the analytic ones obtained for non-random case [6]. 



Figure 1: S^zii^iOj) vs. oj for several values of re. 

In Fig. 1 we plotted the frequency dependence of the random-averaged dynamic 
structure factor for several values of wave vector re = | (a), re = | (b), re = ^ (c). 
One can easily see the differences between the cases a — —1.01 {1), a — 1.01 (2), and 
the case when the transverse field at site is random Lorentzian variable independent 
on the values of exchange couplings at neighbouring sites (3). The frequency profiles 
for non-random case F = are depicted by dashed lines (4). The well- pronounced 
differences in the plotted frequency shapes can be explained bearing in mind that zz 
dynamics is conditioned by excitation of two magnons and taking into consideration 
the changes in the density of magnon states caused by the correlated disorder with 
opposite signs of a and the not correlated disorder. 
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